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Nucler.r  resonance  profiling  was  used  to  measure  the  distributions  of  hydrogen  incorporated  into 
drv  ;;i0„  films  by  thermal  treatments  in  steam.  Thermal  oxides  were  grown  on  silicon  to  a 
tl  ickntss  of  260  nm  in  dry  02  and  were  subsequently  treated  in  steam  at  temperatures  of  320  and 
500C  for  periods  lasting  between  390  and  6  X  105  s.  The  concentrations  of  hydrogen  carried  in  by 
permeating  water  were  then  profiled  with  6.4  MeV  ,SN  ions  using  the  resonant  nuclear  reaction 
'H(I5N,  ay)'2C.  Water  was  seen  to  penetrate  the  films  rapidly  and  to  slowly  react  with  the  Si02 
uniformly  throughout  the  films.  Two  distinct  stages  were  observed  in  the  buildup  of  H,  indicating 
that  the  water/Si02  reaction  involves  at  least  two  concurrent  processes  rather  than  a  single-stage 
process. 

PACS  numbers:  66.30.  —  h,  68.60.  +  q,  81.60.Dq 


Thermal  oxidation  of  silicon  to  form  high-quality  insu- 
ating  films  of  Si02  is  an  important  process  in  modem  elec¬ 
tronic  technology.  The  properties  of  the  oxides  formed  by 
this  process  have  been  found  to  depend  critically  on  the  pres¬ 
ence  of  water  absorbed  in  the  films  either  during  oxidation  or 
in  subsequent  high  temperature  processing.  An  example  of 
this  is  the  recent  observation 1  of  morphological  differences 
between  wet  and  dry  ultra-thin  metal-oxide-semiconductor 
(MOS)  gate  oxides:  wet  films  were  seen  to  be  free  of  the  mi¬ 
cropores  present  in  dry  films  and  showed  better  breakdown 
resistance.  Physical  models  of  the  diffusion  and  incorpora¬ 
tion  of  water  in  SiO,  have  received  comparatively  little  ex¬ 
perimental  testing.  In  an  early  study  of  water  diffusion  in 
fused  silica,  Moulson  and  Roberts2  proposed  that  the  trans¬ 
port  involved  mobile  OH  groups  arising  from  the  reaction 

H20  +  Si-0-Si*=t2SiOH ,  ( 1 ) 

at  the  outer  surface.  They  modeled  the  process  with  Fick’s 
second  equation  using  a  constant  effective  diffusion 
coefficient 

d  [OH]/dt  —  .£>,,* V2 [ OH ] ,  (2) 

determining  Dctt  by  observing  the  time  behavior  of  the  total 
OH  content  during  the  early  stages  of  steam  absorption  and 
vacuum  desorption.  An  alternative  model,  proposed  by  Dor- 
emus,3  was  the  transport  of  dissolved  molecular  water 
through  “doorways”  between  cells  or  interstices  in  the 
network,  with  the  SiOH  groups  produced  by  reaction  (1) 
being  immobile.  He  also  assumed  (a)  that  the  reaction  pro¬ 
ceeded  rapidly  compared  to  diffusion,  allowing  the  equilibri¬ 
um  concentrations  of  diffusing  water  and  SiOH  groups  to  be 
related  by  the  law  of  detailed  balance 

[H,0]  =  [OH]2/*2,  (3) 

and  (b)  that  the  water  concentration  and  its  time  derivative 
were  negligible  compared  to  those  for  SiOH.  These  assump¬ 
tions  led  him  to  model  the  process  with  the  concentration- 
dependent  diffusion  equation 

d  [OH]/*  =  {2D  /K  2)V2[OH]2,  (4) 


where  D  and  Dt„  are  related  by  Dc„  =  2D  [H20]/[OH].  This 
produced  a  much  better  match  to  the  [OH]  profiles  mea¬ 
sured  by  Roberts  and  Roberts4  in  fused  silica  than  is  possible 
with  a  solution  to  Eq.  (2).  (A  profile  measured  by  Moulson 
and  Roberts  was  insufficiently  precise  to  distinguish  the  two 
models.)  This  model,  unlike  the  earlier  one,  also  implies  a 
larger  effective  diffusivity  for  absorption  than  for  desorp¬ 
tion,  consistent  with  the  measurements  of  Moulson  and  Ro¬ 
berts.  Doremus  also  used  this  model  to  treat  the  steady-state 
steam  oxidation  of  silicon,5  relating  the  temperature  and 
pressure  dependence  of  the  parabolic  oxidation-rate  con¬ 
stant  to  that  of  O',, . 

Doremus’  model  applies  to  the  regime  where  equilibri¬ 
um  has  been  established  between  [HzO]  and  [OH].  Our  pre¬ 
vious  observations6  of  tracer  l80  exchange  during  water  dif¬ 
fusion  at  relatively  low  temperatures  were  generally 
consistent  with  this  model.  However,  we  saw  evidence  of  a 
systematic  trend  toward  anomalously  low  exchange  rates 
during  the  early  stages  of  absorption  in  dry  oxides.  This  sug¬ 
gested  that  the  [OH]  had  not  reached  its  saturation  value, 
contrary  to  our  expectation  from  Eq.  (4)  that  it  should  have 
reached  it  on  a  time  scale  of  x\/Dc„,  where  x„  is  the  film 
thickness.  We  therefore  decided  to  investigate  the  early  be¬ 
havior  of  water  absorption  in  dry  oxides. 

Samples  were  prepared  by  forming  thermal  oxide  films 
in  nominally  dry  02  at  1 100  C  on  (100)  surfaces  of  5  (1  cm 
(nominal)  n-type  silicon  to  a  thickness  of  260  nm.  The  sam¬ 
ples  were  subsequently  treated  in  flowing  steam  at  1  atm 
according  to  the  procedure  described  in  Ref.  6.  Formation  of 
additional  oxide  during  the  steam  treatments  was  negligible. 
The  concentrations  of  hydrogen  carried  in  by  absorbed  wa¬ 
ter  were  then  profiled  with  6.4  MeV  l5N+  +  ions7  by  detect¬ 
ing  the  characteristic  4.4-MeV  gamma  rays  from  the  reso¬ 
nant  nuclear  reaction  1  H( 1 5N,  ay) 1 2C.  The  beam  current  was 
about  40  nA  and  the  beam  spot  was  about  3x10  mm,  pro¬ 
ducing  negligible  sample  heating.  Data  collection  lasted 
about  8  min  at  each  beam  energy.  The  integrated  ion-beam 
current  was  24  [tC  at  each  beam  energy. 

Detailed  profiles  of  the  near-surface  regions  showed 
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TABLE  I.  Observed  yield  of  4.4]  gamma  rays  emitted  from  reactions  occurring  at  indicated  depths  in  260-nm  SiO,  films:  number  of  counts  recorded  per  24 
/iC  of  incident  ”N  ions. 


Sample 

Treatment 

temp. 

IC) 

Time 

Is) 

Counts  at  depths 
Inm) 

in  films 

80 

146 

199 

229 

2J8 

Itt 

201 

500 

96  030 

329 

280 

283 

256 

242 

78 

202 

500 

29  160 

256 

211 

221 

219 

208 

43 

203 

500 

14  160 

250 

201 

206 

211 

45 

204 

500 

4  350 

241 

172 

202 

159 

32 

205 

500 

1  860 

209 

193 

139 

114 

43 

301 

320 

601  320 

382 

329 

328 

306 

243 

51 

302 

320 

167  700 

286 

268 

240 

221 

162 

41 

303 

320 

57  600 

217 

235 

212 

164 

147 

43 

304 

320 

21  600 

249 

216 

179 

165 

90 

30 

305 

320 

7  320 

210 

198 

175 

132 

48 

306 

320 

2  760 

218 

193 

151 

94 

46 

307 

320 

1020 

145 

131 

94 

82 

48 

308 

320 

390 

121 

114 

86 

69 

34 

309 

Untreated 

131 

72* 

30 

‘At  depth  of  170  nm. 


that  a  hydrogenous  layer  about  5  nm  thick  coated  the  surface 
of  every  sample,  including  an  unoxidized  Si  sample  which 
was  profiled  to  determine  instrumental  background.  This 
layer,  which  commonly  forms  on  Si02  surfaces  exposed  to 
air,  allowed  us  to  determine  the  effects  of  beam  spread,  ener¬ 
gy  straggling,  and  resonance  width  on  the  depth  resolution 
of  the  profile.  A  yield  curve  of  gamma  rays  collected  from 
the  layer  was  matched  to  a  theoretically  obtained  curve  gen¬ 
erated  by  assuming  a  Lorentzian  resonance  in  the  cross  sec¬ 
tion.  The  best  fit  to  its  effective  width  was  about  8  keV,  corre¬ 
sponding  to  a  depth  resolution  of  about  5  nm  near  the 
surface;  energy  straggling  approximately  doubled  this  near 
the  oxide-silicon  interface.  The  [H]  was  therefore  sampled 
only  every  30  nm  and  no  closer  than  80  nm  to  the  surface  to 
avoid  contributions  from  the  hydrogenous  outer  layer. 

Some  difficulties  have  been  reported  by  other  investiga¬ 
tors  involving  changes  in  H  profiles  induced  by  nuclear  pro¬ 
filing;  these  have  included  hydrogen  effusion8  from  a-Si  and 
room-temperature  vacuum  desorption  of  water  from  some 
soda-lime  glasses.9  To  check  for  these  effects,  data  collection 
from  points  sampled  early  was  repeated  after  completion  of 
profiling.  No  statistically  significant  variations  in  [H]  were 
observed,  either  in  the  interior  or  at  the  surface. 

Our  profiling  results  are  listed  in  Table  I,  which  shows 
the  number  of  counts  recorded  from  reactions  occurring  at 
various  depths  below  the  base  of  the  hydrogenous  layer 
(depths  were  calculated  using  a  specific  energy  loss  of  1.69 
MeV/^m).  The  samples  are  in  two  groups,  one  whose  treat¬ 
ments  were  at  500  C  and  the  other  at  320  C;  also  shown  is  an 
oxidized  but  untreated  sample  (number  309).  An  unoxidized 
silicon  sample  (not  shown)  was  included  for  background  de¬ 
termination;  this  yielded  about  36  counts  at  every  depth. 
After  background  subtraction,  the  counts  were  converted  to 
H  concentrations  with  calibration  data  obtained  from  pre¬ 
vious  measurements  of  samples  containing  known  quantities 
of  hydrogen.10  Some  representative  profiles  are  plotted  in 
Fig.  1;  the  error  b<-'  •’dicate  statistical  uncertainty. 

These  results  show  that  the  [H]  profiles  are  essentially 
uniform  throughout  the  oxides,  rather  than  taking  the  form 
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of  advancing  fronts  of  diffusant.  This  is  clearly  contrary  to 
the  prediction  of  the  Fickian  OH-diffusion  model,  Eq.  (2).  In 
the  context  of  the  molecular  H20  diffusion/reaction  model, 
it  indicates  that  chemical  equilibrium  had  not  been  achieved 
during  the  permeation  process,  invalidating  Eq.  (4)  as  well. 
This  can  be  seen  by  inspection  of  Fig.  2,  in  which  are  plotted 
the  profiles  predicted  by  solving  those  equations.  The  rela¬ 
tive  flatness  of  the  observed  profiles  is  an  indication  that  the 
diffusion  through  the  oxides  was  rapid  compared  to  the  reac¬ 
tion,  allowing  the  diffusant  concentration  to  approach  uni¬ 
formity  before  any  appreciable  quantity  of  product  had  built 
up.  The  one  exception  observed  was  the  oxidized  but  un¬ 
treated  sample,  which  shows  a  small  but  detectable  accumu¬ 
lation  of  H  near  the  surface;  this  may  have  entered  the  oxide 
immediately  after  its  formation. 

The  relative  uniformity  of  the  distributions  under  the 
conditions  of  this  experiment  (i.e.,  thin  films  and  moderate 


FIG.  I.  Representative  [H]  profiles  measured  with  6.4-MeV  "N  ions 
for  some  samples  listed  in  Table  I.  All  samples  were  260-nm  thick  ox¬ 
ides  which  had  been  thermally  grown  in  nominally  dry  Oj  at  1200  C 
and  then  treated  in  ]20  C  steam  at  1  atm  for  the  durations  indicated. 
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FIG  2.  Theoretical  [H]  profile*  at 
various  times  for  diffusion  into  an 
initially  dry  oxide  film  of  thick¬ 
ness  Xa,  normalized  to  [H]  at  infi¬ 
nite  time.  Curve  labels  denote 
time,  which  is  expressed  in  units 
of  x*0/D  The  profiles  shown  in  (a) 
were  obtained  by  solving  Eq.  (2), 
and  those  in  |b)  were  obtained  by 
solving  Eq.  (4)  with  K  =  2;  the 
boundary  conditions  were  in  both 
cases  [H]  (surface)  =  1  and 
3  [H]/3x  (interface)  =  0.  They  may 
be  compared  to  the  corresponding 
profiles  given  by  Crank  for  an  in¬ 
finite  medium  (see  Ref.  11). 


temperatures)  enabled  the  hydration  process  to  be  singled 
out  for  study  independently  of  diffusion.  For  this  purpose, 
the  measured  data  was  summarized  by  averaging  for  each 
steam  treatment  the  [H]  in  the  interiors  of  the  oxide  films. 
Results  of  this  procedure  appear  in  Fig.  3,  which  may  be 
regarded  in  this  context  as  an  example  of  a  Powell-plot12 
used  for  reaction-rate  law  determination.  Note  that  a  break 
clearly  occurs  in  the  curve  of  [H]  vs  t  for  either  treatment 
temperature,  implying  that  the  hydration  of  Si02  involves  at 
least  a  dual  rather  than  a  single  reaction.  One  component 
reaction  appears  to  approach  completion  before  the  other 
has  contributed  significantly:  the  value  of  [H]  arising  from 
the  fast  component  appears  to  saturate  at  about  1.5  X  102° 
cm  ~  3  at  320  C.  This  is  well  below  our  projection  of  4. 1  X  102° 
cm'3  for  the  equilibrium  concentration  at  that  temperature, 
which  we  obtained  by  extrapolating  the  plot  of  OH  solubility 
versus  temperature  calculated  by  Shackelford  and  Ma- 
saryk.13  (We  are  aware  that  other  investigators  have  ob¬ 
tained  different  results.)  At  500  C  the  fast  component  ap¬ 
pears  to  saturate  at  2  X  1020  cm'1,  accounting  for  most  of  the 
projected  equilibrium  concentration  of  2.5 X  1020  cm"’.  A 
rough  estimate  of  the  activation  energies  can  be  obtained 
from  the  buildup  rates  for  the  products  of  the  component 


FIG.  3.  Mean  (H)  observed  in  steam-treated  oxide  films  as  functions  of 
treatment  duration  at  fixed  temperatures  indicated.  (The  curves  con¬ 
necting  the  points  are  intended  only  as  viewing  guides.) 
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reactions  at  different  temperatures.  This  procedure  yielded 
respective  activation  energies  of  approximately  0-2  and  15- 
18  keal/mol  for  the  fast  and  slow  components,  respectively. 

Although  the  details  of  the  hydration  process  are  not 
yet  fully  clear,  the  results  of  this  and  related  experiments 
allow  some  inferences  to  be  drawn.  Our  evidence  for  a  multi- 
component  reaction  is  consistent  with  the  results  of  Walra- 
fen  and  Samanta,14  who  showed  that  the  3690-cm'1  OH 
stretching  contour  which  they  observed  in  the  Raman  and 
infrared  absorption  spectra  of  fused  silica  could  be  decom¬ 
posed  into  as  many  as  four  distinct  components.  The  low 
activation  energy  which  we  infer  for  our  fast  component  is 
consistent  with  their  suggestion  that  hydrogen  bonding 
plays  an  important  role  in  the  reaction;  its  magnitude  is  typi¬ 
cal  of  those  associated  with  that  type  of  bonding. 15  The  com¬ 
parative  rapidity  of  that  component  reaction  and  its  limited 
extent  suggest  the  possibility  that  it  represents  hydrogen 
bonding  between  network  oxygen  and  molecular  water.  The 
higher  activation  energy  of  the  slow  component  suggests  its 
interpretation  as  the  dissociation  reaction  which  Mikkel- 
sen16  recently  observed  to  produce  complete  oxygen  ex¬ 
change  between  network  Si02  and  diffusing  water,  exclud¬ 
ing  the  possibility  of  a  limited  number  of  fixed  saturable  sites 
for  that  reaction.  Further  clarification  of  the  hydration 
mechanism  awaits  more  detailed  studies  of  reaction-product 
evolution.  These  are  now  in  progress,  as  are  model  calcula¬ 
tions  of  total  H  buildup  and  ,sO  exchange  kinetics. 
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